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JOST FUNCTIONS AND JOST SOLUTIONS FOR 
JACOBI MATRICES, II. DECAY AND ANALYTICITY 

DAVID DAMANIK 1 ' 2 AND BARRY SIMON 1 ' 3 

Abstract. We present necessary and sufficient conditions on the 
Jost function for the corresponding Jacobi parameters a n — 1 and 
b n to have a given degree of exponential decay. 
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1. Introduction 

Among the most interesting results in spectral theory are those that 
give equivalent sets of conditions — one set involving recursion coef- 
ficients and the other involving spectral data. Examples are Verblun- 
sky's version [2E] of the Szego theorem (see [22J ) , the strong Szego 
theorem written as a sum rule (see [221 )> the Killip-Simon theorem jTB] 
characterizing L 2 perturbations of the free Jacobi matrix, and Baxter's 
theorem [T1I22]. 

Our goal in this paper is to present such an equivalence for Jacobi 
matrices concerning exponential decay. That is, we consider orthogonal 
polynomials on the real line (OPRL) whose recursion relation is 

xp n (x) = a n+ ip n+1 (x) + b n+1 p n (x) + a n p n -i(x) (1.1) 

for Jacobi parameters {a n }^L 1; {b n }^ =l . Here p n (x) are the orthonor- 
mal polynomials and P-i(x) = (i.e., a is not needed in (jl.lj) for 
n = 0). 

(jl.lj) is often summarized by the Jacobi matrix 
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2 D. DAMANIK AND B. SIMON 

The model of what we will find here is the following result of Nevai- 
Totik i n the theory of orthogonal polynomials on the unit circle 
(OPUC): 

Theorem 1.1 (Nevai-Totik [Ej; see Section 7.1 of Let dfi be a 

probability measure on <9B obeying 

d9 

dfi = w{6) — + dfi s (1.3) 
Fix R > 1 . Then the following are equivalent: 

(1) The Szego condition holds, dfi s = 0, and the Szego function, D(z), 
has D(z)^ 1 analytic in {z \ \z\ < R}. 

(2) 

limsup|a n | 1/n < fT 1 (1.4) 

n^oo 

The Szego condition is 

J log(w(9)) ^ > -oo (1.5) 

in which case D is defined initially on D by 

D(,)=exp^-^-log(^))-j (1.6) 

In (II. 4|) . a n are the Verblunsky coefficients, that is, the recursion coef- 
ficients for the monic OPUC, $ n , 

&n+l(z) = Z$ n (z) - a n $* n (z) (1.7) 

with 

$;(z) = z"$^I7i) (1.8) 

See |22l E31 ESI 021 HH for background on OPUC. 

Also relevant to our motivation is the following simple result: 

Theorem 1.2. Let dfi be a probability measure on <9B obeying (I1.3J) . 
Then the following are equivalent: 

(1) The Szego condition holds, dfi s = 0, and the Szego function, D(z), 
has D(z)~ l a polynomial of exact degree n. 

(2) aj = for j > n and a n ^\ ^ 0. 

Proof. (2) =>■ (1). In this case (see |22[ Theorem 1.7.8]), dfi = 

|y^(e i(9 )|~ 2 , so D^ 1 = y?* is a polynomial. 
(1) => (2). D{z)~ 1 is nonvanishing on D, and so the measure has the 
form §\D(z)\ 2 and so has ctj = for j > n ([221 Theorem 1.7.8]). 
Thus D(z)^ 1 = y?* (z), and since ip n has degree exactly n, $* = — 
ct n _iz$ n _i implies a n _i ^ 0. □ 
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In our work, the spectral measure has the form 

dj(x) = f{x) dx + dj s (1.9) 

where supp/ C [—2,2]. We say d^ s is regular if 7 S ([— 2,2]) = and 
d% has finite support (i.e., no embedded singular spectrum and only 
finitely many bound states). The m-function associated to dj is defined 
on C\supp(c?7) by 

and M is defined on D = {z \ \z\ < 1} by 

M(z) = -m(z + z- 1 ) (1.11) 

Since z t— > z + z~ l maps D to C U {oo}\[— 2,2], M is analytic on 
D\{z G RflD | z + z~ l is a point mass of d'-f} with simple poles at the 
missing points. 

The Jost function, u(z), is defined and analytic on D in many cases 
and determined first by 

\u{e ie )\ 2 lmM(e ie ) = sin9 (1.12) 

where the functions at e %e are a.e. limits as r f 1 of the functions at 
re 10 . The second condition on u is that, for zGi, 

u(z) = z + z^ 1 is a point mass of d'j (1.13) 

If one has the sufficient regularity of Im M on <9D and 7 S is regular, 
(fTTT2l / (ITTT3J) determine u via 

In addition, if the Jacobi parameters obey 



^2 \ a n - 1| + \K\ 



< oo 



n=l 



then the Jost function can be directly constructed using variation of pa- 
rameters (see Teschl |27j). perturbation determinants (see Killip-Simon 
[To]), or an approach of Geronimo-Case 9J. Since this latter approach 
is not well-known and those authors do not provide the detailed esti- 
mates we will need, we have described this approach in Appendix A. 

When there are zeros of u in D, then u does not uniquely determine 
d'y. f is determined by (jl,12j) and 

/(2cos#) = lmM(e id ) (1.15) 
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and the positions of the point masses are the zeros, but the weights, 
Wj, of the zeros (i.e., the values of j({Ej}) = Wj) are needed. The 
possible values of Wj are constrained by 

r sin 2 9 ,„ 

^ Wl + 2 l Wv M = l (116) 

by (fTTHjl . JTT2J), and 

/2 /*7T 
f(E)dE = 2 / /(2cos0)sin0«20 (1.17) 
-2 JO 

Thus, modulo some regularity issues, the knowledge of a c?7 with 
regular djs is equivalent to the knowledge of u and the finite number 
of weights Wj constrained by 1)1.16)1 . Our main goal in this paper is 
to describe what Jost functions and weights are associated to a n 's and 
6 ra 's with a given rate of exponential decay or with finite support. We 
will view the Jost function/weights as spectral data. This is justified 
by the following: 

Theorem 1.3. Let u be a function analytic in a neighborhood of D 
whose only zeros in this neighborhood lie inDflK with those zeros all 
simple. For each zero in B n R, let a weight Wj > be given so that 
()1.16j) holds. Then there is a unique measure d'y for which u is the Jost 
function and Wj the weights. 

Since this is peripheral to the main thrust of this paper, we do not 
give a detailed proof, but note several remarks: 

1. Related issues are discussed in Paper I of this series 

2. One first shows that the M defined by d'y has a meromorphic con- 
tinuation to a neighborhood of D; this is done in Theorem 13.7.1 

of [23]. 

3. The methods we use in Sections El and |3J then show that the a n — 1 
and b n decay exponentially. 

4. Thus, by the results of Appendix A, a Jost function, u, exists. 
u/u has removable singularities, is nonzero on D, is analytic in a 
neighborhood of D, and on <90, \u/u\ = 1. Thus, u = u. 

The perturbation determinant can be defined by 

- U M (Li8) 

This is obviously normalized by 

L(0) = 1 (1.19) 
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which is simpler than Ijl.lfij) . Of course, u(0) can be recovered from 
{wj}f =1 and L(z) by (TTTHJ) and (Onjl . We note that when J - J is 
trace class, we have (see [Tfi] ) 

= det(l + (J - J )[J - + z- 1 )}- 1 ) (1.20) 

Our goal in this paper is to prove four theorems: two in the simple 
case where there is no point spectrum and two in the general case. In 
each pair, one describes finite support perturbations and one, exponen- 
tial decay. We begin with the case of no bound states: 

Theorem 1.4. If a n — 1 and b n = for large n, then L(z) is a 
polynomial. Conversely, any polynomial L(z) which obeys 

(i) L(z) is nonvanishing on D\{±1} 

(ii) If +1 and/or —1 are zeros, they are simple 

(iii) L(0) = 1 

is the perturbation determinant of a unique Jacobi matrix and it obeys 
a n = 1 and b n = for all large n. 

Remark. By Theorem IA.ll there is a precise relation between the 
degree of L and the range of (a n — 1, b n ). 

Theorem 1.5. Let R > 1. If 

lim (K-ll + KI) 1 / 2 "^- 1 (1.21) 

n— >oo 

then L(z) has an analytic continuation to {z \ \z\ < R}. Conversely, if 
L(z) is analytic in {z \ \z\ < R} and obeys (i)-(iii) from Theorem \l.J\ 
then (TT2T1) holds. 

For bound states, things are more complicated. One way of under- 
standing this is the following. Consider a finite support set of Jacobi 
parameters with I bound states. One can change the I weights without 
changing L(z). It is known in that case that changing a single weight 
introduces an explicit (in terms of solutions of the original Jacobi re- 
cursion) correction which decays exponentially (see [21 pp. 64-66] and 
[T31 IT1] ) . That means one expects only one out of the ^-parameter fam- 
ily of Jacobi matrices to have Jacobi parameters with finite support. 
Put differently, there will be many Jacobi matrices with bound states 
but only exponentially decaying Jacobi parameters that have polyno- 
mial Jost functions. So the key is identifying the weights that single 
out finite support. 

Rather than discuss weights, it is more convenient to use residues of 
poles of M. Of course, if Zj G D, Zj + zj 1 = Ej is a point mass in 7, 
then 

Wj = lim (E - Ej)m(E) = {zj 1 - z^zj 1 lim (z - Zj )M{z) (1.22) 
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so the data are equivalent. 

Definition. Let M be the M-function associated to a u and a set of 
weights. Suppose u is analytic in {z \ \z\ < R} for some R > 1 and 
u(zj) = with \zj\ > Rr 1 . We say the weight at Zj is canonical if and 
only if 

Wj = lim (z - Zj)M(z) = -{zj - z^ 1 )[u'{z j ) u(l/E) ] _1 ( L23 ) 

Here are our main theorems on the general case: 

Theorem 1.6. If a n = 1 and b n = for large n, then L(z) is a 
polynomial and all the weights are canonical. Conversely, if L is a 
polynomial obeying 

(l) L(z) is nonvanishing on D\M 

and (ii)-(iii) of Theorem \l-4\ then there is at most one set of Jacobi 
parameters with a n = 1 and b n = for n large that has that L as 
perturbation determinant. Moreover, the weights associated to this set 
are the canonical ones. If these canonical weights lead to Wj > 0, then 
there is a set of Jacobi parameters with a n = 1 and b n = for large n. 

Remark. It is easy to construct polynomial L's which are not 
the perturbation determinant of any finite support Jacobi parame- 
ters, although they are perturbation determinants. For example, if 
L(z ) = L^Zq 1 ) = for some z G (0, 1), (jl.23J) cannot hold. Thus 

L{z) = {l-2z){l-\z) 

is a perturbation determinant but not for a Jacobi matrix of finite 
support. There are also examples where the canonical weights are 
negative. 

Theorem 1.7. Let R > 1. // (11.21)1 holds, then L(z) has an analytic 
continuation to {z | \z\ < R} and the weights of all Zj with \zj\ > i? -1 
are canonical. Conversely, if L(z) is a function analytic in {z | \z\ < 
R} obeying (i')-(iii) of Theorem \l.b\ then (11.21)1 holds if and only if all 
weights for Zj with \zj\ > R^ 1 are canonical. 

These four theorems have a direct part (i.e., going from {a n , & n }^=i to 
L(z)) and an inverse part. The direct parts (except for the importance 
of canonical weights) are well-known. We provide a proof of all but 
the canonical weights in Appendix A. The canonical weight result is 
proven in Section El 

The inverse parts are more subtle — and the main content of this 
paper. The no bound state results appear in Section El and the bound 
state results in SectionEl Our approach is based on the use of coefficient 
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stripping, that is, relating u, M for {a n , & n }^=i to u, M for {a n , b n }^ =1 
where a n = a n+ i, b n = b n+ x- Section [2] will rely on a remarkably simple 
contraction argument, Section El on the fact that coefficient stripping 
only preserves analyticity if weights are canonical. 

One can wonder if one can't at least prove the no bound state results 
by appealing to the Nevai-Totik theory and the Szego mapping (see 
Section 13.1 of [23]) relating OPUC and OPRL. Indeed, we will show in 
Section[2]that our method can be used to prove the inverse part of their 
result. There is a difficulty with blind use of the Szego map, already 
seen by the fact that Jo does not map into Verblunsky coefficients with 
exponential decay (see Example 13.1.3 of [2H])- This can be understood 
by noting that the Jost function, u, for d'j and the Szego function, D, 
for = Sz -1 (g?7) are related by 

*M- = ^ d.24) 
1 — z z 

Thus, D(z)^ 1 is not analytic where u is, unless u(+l) = u(—l) = 0. In 
that case, one can use Nevai-Totik to obtain Theorems 11.41 and II .51 

There are two strategies for dealing with the general case. First 
(and our original proof), one can add extra a's and 6's at the start 
to produce u(+l) = u(—l) = 0. Second, // i— > Sz(/i) is one of four 
maps (see Section 13.2 of |2Sj)- Sz2 maps onto all Jacobi matrices with 
spectrum on [—2,2] and with u(l) ^ ^ u(— 1) and has no division 
factor. Sz 3 and Sz 4 divide by 1 — z and 1 + z and are onto all matrices 
with u{— 1) 7^ and it(l) 7^ 0. In this way, one can always find a \x 
with 7 = Szj(/x) so the /J-function for \i is analytic. 

It should also be possible to prove the inverse results we need us- 
ing the Marchenko equation. That said, we prefer the approach in 
Section EJ 

Surprisingly, the four main results of this paper appear to be new, 
although for Schrodinger operators with Yukawa potentials, there are 
related results in Newton and Chadan-Sabatier Geronimo |S] 
has a paper closely related to our theme here, but he makes an a priori 
hypothesis about M that means his results are not strictly Jacobi- 
parameter hypotheses on one side. So he does not have our results, 
although it is possible that one can modify his methods to prove them. 

An analog of our results on what are Jost functions for Jacobi ma- 
trices of finitely supported Jacobi parameters is the study of the sets 
of allowed resonance positions for half-line Schrodinger operators with 
compactly supported potentials. There is a large literature on this 
question [3 d HZI CHI I2U 123 EOI In particular in [13 d], Korotyaev 
makes some progress in classifying all Jost functions in this case. 
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We announced the results in pE] and some of them have been pre- 
sented in but we note an error in (22]: Theorem 13.7.4 is wrong 
because, when stating existence of a finite-range solution, it fails to 
require u^zj 1 ) ^ and that the canonical weights be positive. 

2. The Case of No Bound States 

Our goal in this section is to prove Theorems ll.4l and ll.5l We suppose 
we have a set of Jacobi parameters {a n , & ra }^=i with Jost function, 
u(z) = u {0) (z), and M-function, M(z) = M^(z). 

Associated to Jacobi parameters {a>k+n,bk+ n }kLi-i we have corre- 
sponding Jost function, u^ n \z), and M-function, M^ n \z). u^ n \z) is 
the solution of a difference equation at where the solution is asymp- 
totic to z n as n —>■ oo. It follows that 

u n (z) = a~ 1 z n u^(z) (2.1) 

obeys (see (jA.4f)|t ) 

a n u n+ i + (b n - (z + z^))^ + a n -iu n ^i = (2.2) 
Moreover (see (lA~39j) and (|A~42|0 . 

M^{z) = Un+1 [ z \ (2.3) 
a n u n (z) 

This leads to the following set of update formulae: 

u ( n+1 \z) = a n+1 z- 1 u {n \z)M {n \z) (2.4) 

M« (zy 1 =z + z~ l - b n+l - a 2 n+l M^ (z) (2.5) 
Since M(z) = (Sq, (z + z^ 1 — J)~ 1 5q), we see 

^L_^l = i + (z) (2.6) 

so that ()2.5|) implies 

' M^(z) \ _1 

~~~J = 1 - W - K + i - 1> 2 + 0(z*) (2.7) 
which means 

log ^ =b n+1 z+((a 2 n+1 -l) + lbl + JJ + Otf) (2.8) 



There is an additional feature we will need. Suppose u(z) is analytic 
in {z | \z\ < R} for some R > 1. Define 



/"(*) = (2-9) 
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for z G Ar = {z \ R 1 < \z\ < R}. Then we claim 

M(z) - M s (z) = [u{z)u\z)]- x {z - z- 1 ) 
To see this, we note (see f)1.12p ) 



(2.10) 



lmM(e ie ) = [u^uie*)]- 1 ^ 



(2.11) 



This is fl2~TU|) for z = e ie , so (|2~TU1) follows by analyticity. 

The strategy of our proof will be to control vS n \ M^ n > inductively 
using (j23D, (jfSD for z G D and (l2~Tn|) outside D. We will then feed 
this control into (|2.8j) to control a n +i — 1 and 6 n +i- 

We want to use the update equations to confirm that vS n > is analytic 
in at least as big a region as u. This will need the assumption that u 
is nonvanishing on D (and will be the key issue to be addressed in the 
next section). 

Theorem 2.1. If u is analytic in {z | \z\ < R} and nonvanishing on 
D\{+1, —1} with at most simple zeros at ±1, then the same is true of 
each vS n ' . 

Proof. By induction, we only need this for By ()2.10|) . M is mero- 
morphic on {z \ \z\ < R} since we can use (12.10}) to define M(z) as a 
meromorphic function in {z | 1 < \z\ < R} and (|2.11|) says the function 
has matching boundary values on \z\ = 1. 
Moreover, (J2.10j) implies 

(i) M(z) has a pole at z^ with 1 < \z^\ < R only if u(z) has a zero 
there and the order of the pole is the same as the order of the 
zero. This is because vP(z) is nonvanishing near z = Zk and M\z) 
is regular near z = Zk, since z^ 1 G D, and we are supposing no 
bound states. 

(ii) If u vanishes at +1 or —1, M has a first-order pole there (for if 
M{z) has a pole at ±1 with real residue, M$(z) has the opposite 
residue, so M — still has a pole). 

Combining this with (J2.4|) . we see u^(z) is analytic in {z \ \z\ < 
R} for any poles of M are cancelled by zeros of u. Moreover, is 
nonvanishing in D for u is nonvanishing on D\{— 1, 1} and M(z)/z is 
nonvanishing since Im M > on DflC+, so Re M > for z G (0, 1) and 
ReM < for z G (—1, 0). It follows is nonvanishing on 1, 1}. 

And (ii) above shows that even if w(±l) is zero, M has a compensating 
pole. □ 

We will also need 

Theorem 2.2. // the Jost function of Jacobi data {a n ,b n }^ =l has 
finitely many zeros in D and the only zeros on 83 are at ±1 and those 
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are simple, then 

K-1I+ K\ ^0 (2.12) 

and 

M {n \z)^z (2.13) 
uniformly on compacts of D. In particular, for each p < 1, 



sup 

\A<p 



(2.14) 



Proof. Since the weight of the spectral measure is given by (jl.lBj) and 
(I1.12|) . the Szego condition holds and so does the quasi-Szego condition 
of [12]. This plus finite spectrum show Y^=\\ a n ~ 1| 2 + \b n \ 2 < oo by 
the work of Killip-Simon [T^]. Thus (|2.12|) holds. 

That implies the corresponding Jacobi matrix converges in norm 
to Jo so the resolvents converge, which implies ()2.13|) . ()2.14|) is a con- 
sequence of MW (z) /z->1 uniformly. □ 

We now combine (|2.4|) and (|2.1(J|) to write the critical update equa- 
tion: 

u {n+l \z) = a n+1 (l - z- 2 ){u^\z))- 1 + a n+1 z- 2 u^(z)Nl(z) (2.15) 
where 

W.M = ^ (2.16) 

SO 

N*(z) = zM (n)i (z) (2.17) 

()2.15|) looks complicated because of the vtS n ^ term. But consider 
expanding all functions in a Laurent series near {z | \z\ = R\} for 
1 < R\ < R. u( n+1 \z) only has nonnegative powers and and 
thus (1 — z~ 2 )(u^)~ l only has nonpositive powers. Thus the first 
term in (|2.15|) compensates for the negative powers in the second term, 
and that is its only purpose. If we project onto positive powers, it 
disappears! 

We thus define P + to be the projection in L 2 (i?i<9B, |^) onto 
{e in9 }^ =1 . Applying P + to (I2~THD . we find 



u 



(n+l) 



(R ie ie )-u {n+1 \0) = a n+1 P4(R l e ie )- 2 [u^\R 1 e ie )-u^ n \Q)]Nl(R l e i 

(2.18) 

where we used the fact that z 2 u^ n \0)Nf l (z) has only negative Laurent 
coefficients. 
Define 

1 /2 

IIMIk= (/l^(^ie ie )-^(0)| 2 ^ (2.19) 
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for functions analytic in a neighborhood of {z | \z\ < R±}. Since P + is 
a projection in L 2 , we obtain 

||M n+1) Hk < a^^lliVK^e^lUIII^HU, (2.20) 

Proof of Theorem \l.J\ It is easy to see that u is a polynomial of exact 
degree k if and only if u is entire and 

lira G (0, oo) (2.21) 

Thus, by ()2.20j) . if w is a polynomial of degree £ and n > £/2, then 
IH^^HU = 0, that is, ift 1 ' is a constant. But then the weight in M is 
the free one, that is, a J+n = 1, bj +n = for j > 0. Of course, L is a 
polynomial if and only if u is. □ 

Remark. This proof and the direction in the appendix allow us to 
relate the degree of the polynomial u to the support of J — Jo- 

()2.20|) also implies 
Proposition 2.3. For 1 < Ri < R, we have 

limsup |||w (ri) |||5/" < i?r 2 (2-22) 

Proof. Note first that 

sup \N s n (R ie i9 )\ = sup\N n (R^e i9 )\ 
e e 

AfW(z) 



< sup 

\z\<RI x 

which goes to 1 by (I2~T1) . Since a n -> 1 by (f2~T2l. 

(n— 1 ^ 1/n 

Ha i+1 ||JV*(V 
J'=0 

so TTM implies (l2~22l . □ 
Proof of Theorem ] 1.51 If /(^) = Xl^Li Ai-?™, then 

oo 

111/111* = £14. w 

n=l 

is monotone in R%, so (|2.22|) implies, by taking R x — > R, 
limsup |||w^|||i^ < R~ 2 . Since the Cauchy integral formula shows 

sup|/(z)-/(0)|<|||/||| 1+£ 
l*l<i 
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we see that, for every 5 > 0, 

sup \u {n \z) - U W(0)| |i2 - 5\- 2n -> (2.23) 
l*l<i 

which in turn, using tt( n )(0) — > 1 (by ()2.13j) ). implies 

1=2 r , , S1 ^ 2 .1, ^ = , , ^ xlo +0(|i2- 5r 2n ) 
7 |wW(e* 9 )| 2 |wM(0)| 2 Vl 1 7 

which implies 

u ( n \0) = 1 + 0{\R - 5\- 2n ) (2.24) 

Thus the difference between the free weight and the weight for 
is 0(\R-5\- 2n ), so 



lim sup ( sup 

By " 



*l<| 



< iT 2 (2.25) 



limsup|6„| 1/n < FT 2 

limsupKa^-^ + i^l 1 /"^/?- 2 

which implies fll.21|) . □ 

That completes what we want to say about OPRL with no bound 
states. As an aside, we show how the ideas of this section provide 
an alternate to the hard (i.e., inverse spectral) side of the Nevai-Totik 
theorem, Theorem 11.11 Their proof is shorter but relies on a magic 
formula (see (2.4.36) of [22]) 

dy* = a n = - Koo J $ n+1 (e w ) D^ )- 1 dfi(9) 

Our proof will exploit or develop the relative Szego function, SqD, of 
Section 2.9 of Our goal is to prove 

Theorem 2.4. Let d/i be a measure on 83 with dfi s = and so that the 
Szego condition holds. Suppose D(z)' 1 has an analytic continuation to 
{z | \z\ < R} for some R > 1. Then 

lim sup \a n \ 1/n < RT 1 (2.26) 

n— >oo 

So we suppose the Szego condition holds, which is equivalent to 

oo 

^|ct„,| 2 <oo (2.27) 

n=0 
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Let dfjL n be the measure with Verblunsky coefficients {oik+n\kLo and 
its Szego function. F^> and f^ 71 ' are defined by 

F {n \z) = J ^-^d^ n \6) (2.28) 

Geronimus' theorem (see [221 ) says that the relation between the /'s is 
given by the Szego algorithm, 

q n + z/("+l)(z) 

l + a n z/(»+ 1 )(z) 
and the equivalent 

Z f ln+1) ( Z )= f n) -\7n^ ( 2 - 31 ) 

1 - a n f^ n >{z) 

In Section 2.9 of [22], the relative Szego function is defined by (p n = 

(i - Kl 2 ) 1/2 ) 

(f)„£>)(2) = () 2.32 



and it is proven that 

D^(z) 

(*nD)(z) = ^J^L (2 .33) 

which we write as 

D {n+l \z)- 1 = D^ n \z)-\5 n D){z) (2.34) 
It will be useful to rewrite (|2.32|) using (|2.31|) to get 

(<5 " D)(2) = iRii (2 - 35) 



Using 



one finds 



1 F^fzl - 1 
f<ri( z ) = ± LW (2.36) 



{5 n D){z) = \z- 1 M^ n \z) (2.37) 

where 

MW(z) = z(l + a n )(F (n) (z) + 1) - (1 + ^(^(z) - 1) (2.38) 

Interestingly enough, fW(2) for n = appears in the theory of mini- 
mal Caratheodory functions on the hyperelliptic Riemann surfaces that 
occur in the analysis of OPUC with periodic Verblunsky coefficients 
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(see (11.7.76) in [2S]); a related function appears in Geronimo- Johnson 
PH. While [22123 introduced both 5 D and M^(z), its author ap- 
pears not to have realized the relation (j2.37|) . 5 n D is nonsingular at 
z = since M^(z) = at z = (since fW(0) = 1). We note that 
where Section 11.7 of [23] uses M(z), it could use (5oD)(z). The differ- 
ence is the at 0+ is moved to oo_ and the pole at oo + to 0_. We note 
that the relation between M and 5qD is hinted at in ()2.32|) . In gaps 
in supp(d/i) in <9IB>, 5 D has poles at zeros of 1 — zf and zeros at zeros 
of 1 — zf\. By (I2.29|) . 5 D has poles at poles of F and zeros at poles 
of F^\ which is the critical property that M needs in the analysis of 

Section 11. 7 of [23- 

We will also need the analytic continuation of 

ReF^\e ld ) = \D ( - n \e i9 )\ 2 (2.39) 

namely, 

i?(n) + (irW)« = 2D( n )p( n ))* (2.40) 
where Jj is given by (|2.9j) . 

Theorem 2.5. Let R > 1. If ' D~ x is analytic in {z \ \z\ < R}, then F 
is meromorphic there with singularities precisely at the zeros of D^ 1 . 
The order of any pole of F is precisely the same as the order of the zero 
of D~ x . (5qD)(z) is meromorphic in the region with poles precisely at 
the poles of F with order no greater than those of F. (D^)" 1 is analytic 
in {z | \z\ < R} and thus, by induction, so is each (D^)^ 1 . 

Proof. By (I2.4(J|) . we can use 

F = -F s + 2DD i (2.41) 

to define a function meromorphic in {z \ R" 1 < z < R} with poles 
in {z | 1 < z < R} precisely at the zeros of D~ l . Since KeF(e ie ) = 
w{e td ) = \D(e w )\ 2 , 

F(e id ) = -F(Wj + 2 Re F(e id ) = F(e id ) 

F and F agree there and so F extends F to a meromorphic function 
in the required region. 

By ()2.37|) and (|2.38j) . (5qD)(z) is meromorphic in the same region 
with poles of order no greater than those of F. Thus D~ 1 (5qD) is 



analytic in {z \ \z\ < R}. □ 
Now define 

A(z) = \ [z(l + a„)(l - F\z)) + (1 + a )(l + F\z))] (2.42) 

B(z) = (l + a )D s - z^il + a^D 6 (2.43) 
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so, by (J2I3ZJ), ([SinHD , and (jOHjl . 

5 L> = z- 1 ^ + 5L> (2.44) 

and thus 

[D {1 Y l = z- l A(D- 1 - D(0)- 1 ) + B + z~ l AD(0) (2.45) 

Proof of Theorem \2.4\ For functions continuous on {z | \z\ = -Ri}, de- 
fine P + to be the projection onto positive Fourier terms and \\\g\\\ = 
||-P+<?||l 2 - Then ()2.45|) implies, for any R± < R, 

Hl^r'llk < ^r 1 sup \A{z)\ Hl^llk (2.46) 

\z\=Ri 

By induction, 

1 1 1 [D^]- 1 WW < Rl n a ■ ■ ■ an-i 1 1 1 D- 1 1 1 \ Rl (2.47) 

where 

a j= sup \A®(z)\ < | (H-loTj-l) sup (\(1-F®(z))z- 1 \+\1+F®(z)\) 
\z\=Ri ki^r 1 

(2.48) 

Since ^^L l a ™l 2 < °°> we nave suP^-cr- 1 \f i z )\ — * 0' "which im- 
plies, by f)2.36p . that 

sup |(1 — F { - j \z))z~ l \ — > 

and thus, by (|2*1%J) . 

limsup |aj| < 1 

so (|2.47|) implies 

limsup Hl^nilX"^^ 1 

71— +00 

By analyticity, this implies (taking Ri — ► i?) 

limsup |K - l||^ n < iT 1 

n— >oo 

which implies 

lim \a n \ 1/n < R' 1 

n— >oo 

since 

„ de 



a n = I e - l °(w n (6)-l)— □ 
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3. The General Case 

In this section, we will prove Theorems 11.61 and 11.71 A key piece of 
the proofs is the following well-known result: 

Theorem 3.1. If {a n ,b n }™ =1 are a set of Jacobi parameters with 
spec( J)\[— 2, 2] finite, then there is a k so that {a n+ k,b n+ k\ c £ =l are 
a set of Jacobi parameters with spec(Jfc) C [—2, 2]. 

Proof. By a Sturm oscillation theorem (see [2U E]), the number of 
spectral points in (2, oo) is the number of sign flips of {P n (z)}™ =1 and 
in (— oo, — 2) of {(— l) n P n (— 2)}^. By assumption, these numbers 
are finite, so for some k, {P n+fc -i(2)}£° =1 and {(-l) ri P n+fc _i(-2)}~ =1 
have fixed signs. By a comparison theorem and the oscillation theorem 
again, it follows that spec(Jfc) C [—2,2]. □ 

The reason weights have to be canonical is that Theorem 12 . II can fail 
if there are zeros in D. This is because vS 1 ' = a\Z~ l uM and 

M{z) = M\z) + (z- ar 1 )^)^)] -1 (3.1) 

If Zj is a (real) zero in D with \zj\ > R~ l , M(z) may have a pole at 1/ Zj 
due to the (w') -1 term. uM will then have a pole (unless u has a zero, 
which we will see does not help). The way to avoid this is to arrange 
for M* to have a compensating pole, and this will happen precisely if 
the weight at Zj is canonical! Here is the detailed result: 

Theorem 3.2. Letu be a Jost function and be analytic in {z \ \z\ < R} 

and suppose there are finitely many zeros {Zj}j =1 in D. Then M is 
meromorphic in {z \ \z\ < R} and the only possible poles of M in 
{z | 1 < \z\ < R} are at points zj 1 where \zj\ > R" 1 . Moreover, 

(i) If u has a zero of order k > 1 at zj 1 , then M has a pole of order 
k + 1 there. 

(ii) If u(z~ 1 ) 0, then M has a pole at zj 1 if and only if the weight 
at Zj is not canonical. 

In particular, = uM is analytic in {z \ \z\ < R} if and only if all 
weights at those Zj with \zj\ > R^ 1 are canonical. 

Proof. The zeros at Zj are simple, so if u has a zero of order k > 1 at 
zj 1 , then (u(z)u^(z))^ 1 has a pole of order k + 1. Since M"(z) has only 
simple poles in {z \ \z\ > 1}, M(z) has a pole of order k + 1. This 
proves (i). 

If u(zj l ) ^ 0, both (z - z'^lu^u^z)}- 1 and M^(z) have simple 
poles at zj 1 . Their residues cancel if and only if (jl.23J) holds. □ 

This gets us one step if the weights are canonical. We get beyond 
that because automatically weights after that are canonical! 
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Theorem 3.3. Let u be a Jost function and be analytic in {z | \z\ < 
R}. Let u {1) = a 1 z~ 1 uM and let M (1) obey (j2~5|l . Then zj E D with 
\zj\ > R^ 1 implies that u^^zj 1 ) ^ and the weight of is canon- 
ical. 

Proof. Zeros of u in D are cancelled by poles of M, so u^ l > has 
zeros precisely at points Zj where M(zj) = 0. Since [M(zj) — 
M»(%)]u«(%)m(%) = ~ Zj - zj 1 ^ 0, M»(%) ^ 0. By (US]), M(5~ 1 ) ^ 
implies M^(z) is regular at zj" . By a small calculation, t*W,M^ 
obey (|2.10|) . so the weight must be canonical. □ 

Remark. There is a potentially puzzling feature of Theorem 13.31 
If stripping a Jacobi parameter pair cannot produce noncanonical 
weights, how can they occur? After all, we can add a parameter pair 
before J and then remove it. The resolution is that adding a param- 
eter pair also shrinks the region of analyticity if J has a noncanonical 
weight. Essentially, noncanonical weights produce poles in a meromor- 
phic with residues which are a "resonance eieenfunction." Lack 

of a canonical weight in is a sign that this resonance eigenfunction 
function vanishes at n, but then it will not at n + 1 or n — 1. 

Proof of Theorem ] 1.61 If some weight is not canonical, is not entire, 
and so the Jacobi parameters cannot have finite support. If all weights 
are positive, one can normalize u so (jl.l6j) holds with the canonical 
weights and obtain a positive weight since the necessary canonical 
weights are all positive. With canonical weights, uM is entire and 
by (jSU) and M(z) = 0(z) at z = 0, we see that M(z) = 0(l/z) at 
z = oo so long as u(z) = 0(z l ) with £ > 2. Thus is a polynomial 
of degree at most 1 less than u. Iterating and using Theorem 13.11 we 
eventually get a polynomial of vP*'(z) with no zeros in D and so, by 
Theorem II .41 a finite-range set of {a n ,b n }. □ 

Proof of Theorem \l.l\ If weights are canonical, we can iterate to a vS k ' 
nonvanishing in D, and use Theorem 11.51 If some weight is not canon- 
ical, is not analytic in {z \ \z\ < R}, so (|1.21|) fails on account of 
Theorem IA.31 □ 

Appendix A. The Geronimo-Case Equations 

In this appendix, we provide a proof of basic facts about the Jost 
functions in the case of regular bounds on the Jacobi coefficients. We 
do this primarily because we want to make propaganda for a lovely set 
of equations of Geronimo-Case [5] which have not yet gotten the at- 
tention they deserve. It is also useful to keep this paper self-contained. 
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More well-known approaches to Jost functions involve finding the Jost 
solution by using variation of parameters about the free solutions (see, 
e.g., Teschl [27J) or, going back to Jost-Pais [T3], as perturbation de- 
terminants (see, e.g., Killip-Simon [To]). 

There is some overlap in our presentation and that of Geronimo- 
Case and Geronimo jSj, but we feel it might be useful to present 
the detailed estimates concisely. And we wish to emphasize the a pri- 
ori derivation of the GC equations (rather than presenting them and 
proving they have asymptotic properties) by identifying what their ip n 
(= z~ n g n {z) for our g n below) is. In any event, we have put this mate- 
rial in an appendix since we regard it as a review. 

Given a Jacobi matrix, as in (jl.2|) . we define Oq = 1 and look at 



solutions (/ n )£Lo °f 







n 



1,2, 



(a.i; 



(A.2) 



a n fn+i + (K — E)f n + a n ^if n _ 

Of course, one solution of this is 

/„ = p n -i(E, J) 

with p n the orthonormal polynomials associated to J. 

If /„ and k n are two sequences, we define their Wronskian (we will 
often drop the J), 

W(f, k; J)(n) = a n (f n+1 k n - f n k n+1 ) (A.3) 

If /, k both solve (jA.ljl . the Wronskian is constant. 

Given a Jacobi matrix J, define Je by setting a e+ i = a i+2 = ■ • ■ = 1 
and bi + i = bi +2 = • • • = (in [ 16 j and j2Hj, Ji is defined also setting 
an — 1; it is different from Ji), that is, 

/h ai o \ 



Je 



a i 




a i 
a 2 



a 2 



a e -i be ai 
a t 1 
10 



u n (z 
z + z~ 



V 

; Ji) solves ()A.1|) for Jt with u n {z\ J £ 



(A.4) 



z' 



) 



if n > i + 1 where 



E. u n is called the Jost solution. The Jost function is 

u(z; J e ) = W(p.^(z + z- 1 ; J t ),u.(z, J t )) (A.5) 
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= u (z; J t ) 

by taking the Wronskian at n = 1. Since ug{z\ 3%) = a^z 1 , taking the 
Wronskian at n = £, we find 

u(z; Ji) = ai(p e (z + z~ x \ J t ) — z l - p t -\{z + z\ J £ )z e+1 ) (A.6) 

at 

Since the parameters of Je and J agree for a m , b m , m = 1, 2, . . . , £, 
we have p m (z, J) = p m (z, Je) for m < £. Thus we define 

g n (z) = z n (p n (z + ^ J - a n zp n ^! \z + J J J (A.7) 

which is a polynomial in z of degree at most z 2n . g n {z) is the Jost 
function for J n , that is, 

g n (z) = u(z] J n ) (A.8) 

It is also natural to define 

c n (z)=z n pjz + ~^j (A.9) 

a polynomial in z of exact degree 2n. 
Clearly, (jA.7|) becomes 

9n+i(z) = c n+ i(z) - a n+1 z 2 c n (z) (A. 10) 

The fundamental equation 

a n+1 p n+1 fz + = f Z + ^ - & n +l^Pn ^ + _ a nVn-\ (z + J 

multiplied by becomes (using (jA.10|) for n — > n — 1) 

a n+1 c n+1 = (z 2 - b n+1 z)c n + g n (A. 11) 

Multiplying ()A.10|) by a n+1 and using ()A.11|) . we find 

a n +x9n+i = [(1 - o?n+i) z2 ~ K+lA C n + 9n (A. 12) 

The last two equations, which we call the Geronimo-Case equations, 
are the fundamental recursion relations with initial conditions 

g (z) = c (z) = 1 (A.13) 

(One checks g = 1 by looking at (jA.ll|) using c = 1 and c\ — -jj-(z + 
- — fox)-) K is natural to define the Geronimo-Case update matrix 
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so fPT]) / (fP2|) become 



r " +A ' U n J Cn ) (A.15) 



9n+l J a n+l \9n 



or, with 



T n = U n U n . 1 ...U 1 (A.16) 

^)= (ai ... „)-XwQ) (A.17) 

If we use unnormalized functions 

C n (;z) = a x ■ ■ ■ a n c n (z) G n (z) = ai . . . a n g n (z) 

we have 

^ T B H (A.18) 



Since we have the normalization conditions C n (z) = z 2n + lower order, 
G n (0) = 1, this has some similarity to the equations for orthogonal 
polynomials on the unit circle |2*2*1 I2HJ 1213 El E], a motivation for 
Geronimo-Case j^j, but there are numerous differences; for example, 
C*(z) = G n (^) for the reversal operations of the Szego recursion but 
not here. 

One consequence of this is immediate since u — g n if ag — 1, bg = 
for £ > n + 1 : 

Theorem A.l. If J — J is finite range, then u(z; J) is a polynomial. 
If ae — 1, bg — for £ > n + 1 and a n 7^ 1, then deg(w) = In; if a n — 1 
6u£ 6 n 7^ ; deg(-u) = 2n — 1. 

Remark. Indeed, in the case a n 7^ 1, the proof shows that 

u(« ; J) = ^ ~ a ^ z 2n + lower order (A. 19) 

a x . . . a n 

and if a n = 1 but 6 n 7^ 0, 

u(*; J) = — z 2n ~ l + lower order (A.20) 

ai . . . a n 

Proof. We have, if at = 1, bi = for £ > n + 1, that 

u(z; J) = ^[(l - a n)^ 2 ~ 6n*]cn-i(z) + <fo-i(*) 

Since deg(g n _i) < 2n — 2 and c n _i(z) = (ai . . . a n _ 1 )~ 1 z 2n_2 + lower 
order, the result is immediate. □ 
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Notice if b n = 0, a n = 1, then 

U n = UM(z)=(* }) (A.21) 

We will make one of three successively stronger hypotheses on the 
Jacobi coefficients: 

oo 

^[|& n | + |e4-l|]<oo (A.22) 

n=l 
oo 

Y,n[\b n \ + \a 2 n -1|] <oo (A.23) 

n=l 

\b n \ + \a 2 n - 1| < CR~ 2n (A.24) 

for some R > 1. a 2 — 1 will enter in estimates, but since |a„ — 1| < 
|a^ — 1| < (1 + sup n |a n |)(|a n — 1|) in all these estimates, \a n — 1| can 
replace \a% l — 1| with no change. 

In all cases, J^|a n — 1| < oo, so YYj=i a j * s uniformly bounded above 
and below, and hence c (resp., g) and C (resp., G) are comparable. 

We will first prove bounds and use them to control convergence: 

Theorem A.2. (i) Let (fX"22"jl hold. Then for each z <G B\{±1}, 

sup[|G n (z)| + \C n {z)\) = A (z) < oo (A.25) 

n 

where Aq{z) is bounded uniformly on compact subsets o/D\{±l}. 

(ii) Let (jA.23|) hold. Then for some constant A\, 

sup_\G n (z)\ < A ± (A.26) 

n,z£B 
n,zgB 

(iii) Let f)A.24|) hold and let R > 1. Then there is some constant A 2 
such that for all z with \z\ < R, 

\G n (z)\ + \C n (z)\ < A 2 (l+n)[max(l,|z|)] 2m (A.28) 
Proof, (i) If z ^ ±1, UW(z) can be diagonalized; explicitly, 

L{z)U®\z)L{z)-*= (A.29) 

where 



(A.30) 
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Thus for 2 ^ ±1, 

L{z)T n (z)L(z)^ 1 = (K (z) + B n (z)) . . . (K (z) + B{z)) 

with 



Ko(z) 



z 2 
1 



B n (z) = L{z) 



-b n z 
-b n z + (1 - a 2 n )z 2 



Liz)- 1 



If \z\ <l,z^ ±1, 



\\K (z) + B n (z)\\ < 1 + \\L(z)\\ \\L(z)-i\\ [|1 - al\ + I* 



so (lA~22j) implies ((A~25l) with 

2 OO r 



4o(z) < H 



1 



n 

i=i 



1+ 14 



1 -2 2 



(A.31) 



(ii) Let us show inductively that for all \z\ < 1, 

n 

|G„W|<n[l+j(|l-a3l + N)] 



i=i 



\C n (z)\<(n+l)H[l+j(\l-a 2 \ + \b, 
This implies (|A36|l / (jA~2T|) with 



(A.32) 
(A.33) 



To prove (jA.32j) / ()A.33|) , note they hold for n = (where they say 
I Go | < 1, \C \ < 1). If they hold for n, then by (lA~T2i 

n 

\G n+1 (z)\ < [(n + - 1| + |6 n+1 |) + 1] A [1 + J(|l - a? I + N)] 

and by (1A~TT1) . 

n 

\C n+1 (z)\ < [(n + 1)(1 + + 1] J] [1 +idl " 41 + 

i=i 

(iii) Since U n (z) is analytic, the maximum principle says we need 
only prove the estimate for \z\ = R — e with e small. As in part (i), 

\\K (z) + B n (z)\\ < \z\ 2 [l + \\L(z)\\ Wmi^iKl + |1 - al\)\ 

from which ()A.28|) follows. □ 
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Theorem A.3. (i) Let |X"22|) hold. Then 

u(z; J) = lim g n (z) (A.34) 

converges for all z G D\{±1} uniformly on compact subsets of 
D\{±1}. u is analytic on D, continuous on D\{±1}. Moreover, 
for z GD, we have Szego asymptotics 

lim c n (z) = ^4 (A.35) 

n-+oo 1 — Z 

(ii) Let (jA.23|) hold. Then (jA.34|) holds for all z G D converging 
uniformly there, u is continuous on D. 

(iii) Le£ flA~24j) /iota. T/ien flA~34|) /iotas /or 2 G {z \ \z\ < R] uni- 
formly on compact subsets. 

Remark. In [5 , we proved in general that if (jA.34|) holds, then so 
does (|A~35l) for \z\ < 1. 

Proof. Again, since YYj=i a j ls uniformly bounded above and below, it 
suffices to prove convergence of G n and (jA.35J) in the sense that 

lim CJz) = (1 - z 2 )- 1 lim GJz) (A.36) 

n^oo n— >oo 

By (EH, 

|G n+ i(^) - G„(*)l < (\a 2 n+1 - 1| |^| + |6 n+1 | \z\)\C n (z)\ (A.37) 

Thus in all cases, the previous theorem implies 

00 

^2 \G n +i{z) ~ G n {z)\ < 00 

n=0 

uniformly on compacts of the appropriate regions. 
All that remains is (jA.36|) . We have for \z\ < 1, 

\C n — G n -\ — z 2 C n -\\ < Ao(z)\b n \ 

Iterating we see that 



n—l 



n-1 



Cn-J2 Gn-j-lZ 2 ' < \z\ 2n + A (Z) \ b n 

3=0 



-J I \ Z 



V 



3=0 

For k < n — 1, we have 

n— 1 k—l n—l 

i=o i=o £=fe 
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If we let first n — > oo and then k — > oo, we see that 

|2j 



n-1 
3=0 



n-3 I 







since YlT=o\^ k \ < 00 • Thus 



lim 

71— >OC 



n-1 



C,, 



3=0 



2j 



But we claim 

n-l 

by an inequality like (|A.38|) . Thus 



2j 



1 - z 



2n 







proving ()A.36|) . □ 

To understand where ()A.36j) comes from, the above actually shows 
Too = lim^oo T n exists if z 6 D. Since 

T n+1 = [C/(°)(z) + (CW) - U^\z))]T n 

we see that U<® 0)^0) = T^z) so U ( P{z) k T 00 (z) = r^z). Since 
£A°)(z) has eigenvalues z 2 and 1 and \z 2 \ < 1, we see that 

-P(z)T' 00 Tqq 

with P{z) the spectral projection for eigenvalue 1, that is, TqoQ) is a 
multiple of the eigenvectors of U^°'(z) with eigenvalue 1. This eigen- 
vector is f 1 ^ 1 ] -2 ■*) , explaining why (|A.36J) holds. 

If J is a Jacobi matrix, is the Jacobi matrix with af$ = a n+ £, 
b$ = b n+E . If J obeys any of (jA~22|l - (|A~24|l . so do for all E, and 
thus u(z; J^) exists in the approximate region. We define the Jost 
solution (we will show below it agrees with our earlier definition when 
J = J t ) by 

u n (z; J) = a- l z n u(z; J (n) ) (A.39) 
By the "appropriate region" we mean 

D\{±1} if (lA~22l) holds 

D if (lA~23l) holds 

{z | \z\ < R} if (lA~24j) holds 
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Theorem A. 4. (i) u n (z; J) defined in the appropriate region obeys 
a n u n+1 (z; J) + (b n - (z + z~ l ))u n (z] J) + a n -xU n - X (z\ J) = (A.40) 

(ii) In the appropriate region, 

lim z- n u n (z; J) = 1 (A.41) 

(iii) For z e D, 

u(z; J (1) ) = ai z- x u{z\ J)M(z; J) (A.42) 
where for zGD, 

M(» = + J)- 1 ^) (A.43) 

(iv) T/ie on ft/ zeros of u{z; J) in D are at real points (3 with [3 + (3~ l a 
discrete eigenvalue of J. Each zero of u(z; J) in D is simple. 

(v) The only zeros of u{z\ J) in <9B are possible ones at ±1. If (jA.23|) 
holds, these zeros are simple in the sense that 

u(±l, J) = lim fl-We* J) ^ (A.44) 

(vi) In case (|A.24J) no/ds and m(-2o; J) = u>z£n -R -1 < |zq| < 1> ^en 

(vii) // (|A.22|) holds, M(z) has a continuation from D to D\{±1} ; 
which is everywhere finite and nonzero on <9IB)\{±1} and 

\u{e ie )\ 2 Im M(e ie ) = sin 9 (A.45) 

Proof, (i) Since J^) = (^(z; J) — > u(z; J), it suffices to prove (jA.40|) 
for J = Jg for the new definition of u. Let u n (z; Jg) temporarily denote 
the old definition, that is, the solution of (jA.40|) with u n (z] Jg) = z n 
for large n. If we prove u n (z; Jg) = u n (z; Jg), clearly (|A.40|) holds for 
the new definition. 

Since shifts by k steps and z n = z~ k (z n+k ), we have for all n > 1 
and k > 1, 

u n (z; [Jt\ {k) ) = z~ k u n+k (z; Jg) 

Because uq(z\ [Jtft*') is computed using a = 1 but u n (z; Jg) with a = 
Ofc, we have 

u (z; [Jg] {k) ) = a k z~ k u k (z; Jg) 

Thus 

u k (z; Jg) = a~ l z k u(z; [Jg] {k) ) = u(z; Jg) 
(ii) Clearly, (jA.41|) is equivalent to 

lim Goo(z;J (n) ) = l 
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It follows from |A~37|) that 



oo 

|Goc(^;J (n) )-l|< [^-a^ + fiiMlCj-n&J^l 

j=n+l 

By Theorem I A. 2\ the right-hand side goes to zero as n — > oo. 

(iii) For z G D, u n (z; J) is the Weyl solution of (jA.40jl . that is, the 
unique up to a multiple solution of (|A.40|) which is £ 2 at infinity (since 
(1A.41J) implies u n ~ z n G £ 2 ). Thus, by general principles, the matrix 
elements of R(z) = (z + — J) -1 are 

„/ \ Pmin(n,m) — 1 (^)^max(n,m) \Z) 

W(p, U) 

where W(p, u) is the Wronskian of p, u. Evaluating the W at n = 0, 

W(p, u) = u (z; J) = u(z; J) 
Thus (the Weyl formula for the Weyl-Titchmarsh function) 

R n (z) = ^1 (A.46) 
u (z; J) 

= a " ZU ^ { (A.47) 
u(z; J) 

which is ()A.42|) . 

(iv) Since z~ n u n {z\ J) — > 1 as n — > oo, Mo(^; «/) and «i(z; J) cannot 
have simultaneous zeros. Thus, by (fA~47j) . zeros of m in D are precisely 
poles of M(z), and the order of the zeros is the order of the pole. Since 
M has only simple poles in D and precisely at points with z — /3, 
f3 + (3~ l an eigenvalue of J, u has only simple zeros at the prescribed 
points. 

(v) If z G dD, u n (z; J) and w n (z -1 ; J) solve the same Jacobi equation, 
and so their Wronskian is constant. Since a n — > 1 and u n (z; J) ~ z n , 
the Wronskian at infinity is z — z~ x ] while at zero, it is given in terms 
of Mo and Ui, so 

Ui(z; J^oiz' 1 ; J) - J)u (z; J) = z - z' 1 (A.48) 

for z G <9B. On the other hand, u n (z; J) is real if z is real so, by 
analyticity and continuity, u n (z; J) = u n (z;J). Since z~ l = z for 
zeD, (TQ8|) becomes 

Im(ui(e i0 , J) Uo(e ie } J) ) = sinfl (A.49) 
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implies for 8^0, u {e ie , J) + and if u (±l) = 0, then with 
u' (±l,J) = \im 0lo u o (e id ,J)/(t6), 

- Re(iti(±l, J)u' (±l, J)) = 1 

proving that «(, 7^ 0. 

(vi) u n {z~ x \ J) and w„(z; J) obey the same equation, so (jA.48|) holds. 
This implies it cannot happen that u (z; J) and Uq(z~ 1 ; J) are both 
zero. 

(vii) By (TO! . 

lmM(e w ) = \u (e ie , J)\~ 2 lm( Ul (e ie , J) M (e ie , J) ) 
which, given (jA.49|) . implies ()A.45|) . □ 

Warning. In j^j, we constructed a Jost solution u n {z\ J) for certain 
J's and z G D which may not obey (jA.22|) . It can happen that for such 
J's that u(z; J) has a zero boundary value as z = re td — *> e . This is 
because while z~ n u n (z; J) —>■ 1 for z G D, this may not be true for the 
boundary values. 



References 

[1] G. Baxter, A convergence equivalence related to polynomials orthogonal on 
the unit circle, Trans. Amer. Math. Soc. 99 (1961), 471-487. 

[2] K. Chadan and P.C. Sabatier, Inverse Problems in Quantum Scattering The- 
ory, 2nd ed., Springer- Verlag, New York, 1989. 

[3] E.A. Coddington and N. Levinson, Theory of Ordinary Differential Equations, 
McGraw-Hill, New York- Toronto-London, 1955. 

[4] D. Damanik, R. Killip, and B. Simon, Necessary and sufficient conditions in 
the spectral theory of Jacobi matrices and Schrddinger operators, Int. Math. 
Res. Not. (2004), no. 22, 1087-1097. 

[5] D. Damanik and B. Simon, Jost functions and Jost solutions for Jacobi ma- 
trices, I. A necessary and sufficient condition for Szego asymptotics, preprint. 

[6] L.D. Faddcyev, The inverse problem in the quantum theory of scattering, J. 
Math. Phys. 4 (1963), 72-104. 

[7] R. Froese, Asymptotic distribution of resonances in one dimension, J. Differ- 
ential Equations 137 (1997), 251-272. 

[8] J. S. Geronimo, Scattering theory, orthogonal polynomials, and q-series, SIAM 
J. Math. Anal. 25 (1994), 392-419. 

[9] J.S. Geronimo and K.M. Case, Scattering theory and polynomials orthogonal 
on the real line, Trans. Amer. Math. Soc. 258 (1980), 467-494. 
[10] J.S. Geronimo and R. Johnson, An inverse problem associated with polynomi- 
als orthogonal on the unit circle, Comm. Math. Phys. 193 (1998), 125-150. 
[11] Ya. L. Gcronimus, Polynomials Orthogonal on a Circle and Their Applica- 
tions, Amer. Math. Soc. Translation 1954 (1954), no. 104, 79 pp. 



28 D. DAMANIK AND B. SIMON 

[12] Ya. L. Geronimus, Orthogonal Polynomials: Estimates, Asymptotic Formulas, 

and Series of Polynomials Orthogonal on the Unit Circle and on an Interval, 

Consultants Bureau, New York, 1961. 
[13] F. Gesztesy, A complete spectral characterization of the double commutation 

method, J. Funct. Anal. 117 (1993), 401-446. 
[14] F. Gesztesy and G. Teschl, Commutation methods for Jacobi operators, J. 

Differential Equations 128 (1996), 252-299. 
[15] R. Jost and A. Pais, On the scattering of a particle by a static potential, Phys. 

Rev. 82 (1951), 840-851. 
[16] R. Killip and B. Simon, Sum rules for Jacobi matrices and their applications 

to spectral theory, Ann. of Math. (2) 158 (2003), 253-321. 
[17] E. Korotyaev, Inverse resonance scattering on the half line, Asymptot. Anal. 

37 (2004), 215-226. 

[18] E. Korotyaev, Stability for inverse resonance problem, Int. Math. Res. Not. 
(2004), no. 73, 3927-3936. 

[19] P. Nevai and V. Totik, Orthogonal polynomials and their zeros, Acta Sci. 
Math. (Szeged) 53 (1989), 99-104. 

[20] R.G. Newton, Scattering Theory of Waves and Particles, 2nd ed., Dover Pub- 
lications, Minneola, NY., 2002. 

[21] B. Simon, Resonances in one dimension and Fredholm determinants, J. Funct. 
Anal. 178 (2000), 396-420. 

[22] B. Simon, Orthogonal Polynomials on the Unit Circle, Part 1: Classical The- 
ory, AMS Colloquium Scries, American Mathematical Society, Providence, 
R.I., 2005. 

[23] B. Simon, Orthogonal Polynomials on the Unit Circle, Part 2: Spectral The- 
ory, AMS Colloquium Series, American Mathematical Society, Providence, 
R.I., 2005. 

[24] B. Simon, Sturm oscillation and comparison theorems, Proc. Sturm 200th 

Birthday Conference, Geneva, 2003, preprint. 
[25] B. Simon and A. Zlatos, Sum rules and the Szego condition for orthogonal 

polynomials on the real line, Comm. Math. Phys. 242 (2003), 393-423. 
[26] G. Szego, Orthogonal Polynomials, Amer. Math. Soc. Colloq. Publ., Vol. 23, 

American Mathematical Society, Providence, R.I., 1939; 3rd edition, 1967. 
[27] G. Teschl, Jacobi Operators and Completely Integrable Nonlinear Lattices, 

Mathematical Surveys and Monographs, 72, American Mathematical Society, 

Providence, R.I., 2000. 
[28] S. Verblunsky, On positive harmonic functions (second paper), Proc. London 

Math. Soc. (2) 40 (1936), 290-320. 
[29] M. Zworski, Distribution of poles for scattering on the real line, J. Funct. 

Anal. 73 (1987), 277-296. 
[30] M. Zworski, A remark on isopolar potentials, SIAM J. Math. Anal. 32 (2001), 

1324-1326. 



